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ABKEGE DESCRIPTIF 



Precede de communications radiomobiles AMRT 
iteratif . 

Selon 1' invention, on repartit plusieurs sequences 
de reference dans chaque tranche temporelle allouee. A 
la reception, on a recours a une representation en 
treillis du canal multitra j ets et a un precede 
iteratif . 

Application aux radiocommunicat ions avec les 
mobiles . 

Fig. 12 
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PROCEDE DE COMMUNICATIONS RADIOMOBILES AMRT ITERATIF 
DESCRIPTION 

5 Domaine technique 

La presents invention a pour objet un procede de 
communications radiomobiles AMRT ("acces multiple a 
repartition en temps", en anglais TDMA pour "Time 
Division Multiple Access") iteratif. 
10 Elle peut s'appliquer au systems radiomobile 

europeen GSM [1] et son evolution GSM++ ou au systems 
americain D-AMPS . 

Etat de la technique anterieure 

15 Lss systsmes AMRT permettent de repartir dans le 

temps plusieurs utilisateurs communiquant a travers un 
meme canal radioelect rique [1, 3] . On separe chaque 
utilisateur des autres en lui assignant periodiquement 
une tranche temporelle (ou intervalle, ou creneau 

20 temporel, ou "timeslot" en anglais) pour lui permettre 
d'emettre et une autre tranche temporelle pour lui 
permettre de recevoir. Pour se premunir des 
imperfections en terms de synchronisation temporslls 
entre emstteurs et de 1 ' interference entre tranches 

25 temporelles consecutives produite par la pluralite de 
trajets, des intervalles de garde sont prevus de part 
et d' autre de chaque tranche temporelle. 

La figure 1 annexes montre ainsi des tranches 
temporelles TAi, TA2, TAn allouees a N 

30 utilisateurs, avec des intervallss ds gards IG. 
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Le signal radioelectrique emis durant une tranche 
temporelle est obtenu par transposition en frequence du 
signal en bande de base equivalent. Le signal en bande 
de base est le resultat du filtrage du train de donnees 
5 a emettre par un filtre dit d' emission. Habituellement , 
le train de donnees emis est compose de deux sequences 
de symboles de donnees separees en temps par une 
sequence de symboles connue du recepteur, dite sequence 
de reference. Les deux sequences de symboles de donnees 

10 peuvent eventuellement provenir du codage et de 
1' entrelacement des informations a emettre. 

Afin de permettre une estimation precise du canal 
au niveau du recepteur, la sequence de reference 
possede souvent de bonnes proprietes de correlation 

15 temporelle. A titre d'exemple, les sequences de 
reference utilisees dans le cadre du systerae GSM, 
appelees sequences CAZAC (contraction de "Constant 
Amplitude - Zero Autocorrelation" en anglais) [1, 4, 
5] , sont des sequences a composantes prises dans 

20 1' alphabet bipolaire {-1, 1} et possedant une fonction 
d' autocorrelation circulaire nulle partout sauf a 
1' origine . 

Le canal radiomobile emprunte lors d' une 
communication entre un emetteur et un recepteur est 

25 generalement de type multi-tra j ets avec evanouissements 
rapides dits de RAYLEIGH. L' existence de plusieurs 
trajets vient de ce que I'onde radioelectrique se 
propage selon plusieurs chemins entre le lieu 
d' emission et le lieu de rece^ption. 

30 Le signal re?u est alors la somme de plusieurs 

repliques plus ou moins retardees et plus ou moins 
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alterees en phase et en amplitude. Pour obtenir- une 
restitution fiable de 1 ' information transportee par ces 
repliques, le recepteur opere un filtrage adapte au 
filtre d' emission et au canal et il combine de maniere 
optimale les contributions energetiques de toutes les 
repliques du signal emis. Le signal en sortie du filtre 
adapte est echantillonne au rythme des symboles et il 
est blanchi au moyen d'un filtre discret dit 
blanchisseur . 

Les echantillons en sortie du filtre blanchisseur 
fournissent une version filtree et bruitee du train de 
donnees emis. Le filtre discret associe a cette version 
filtree, egalement appele canal discret, presente une 
reponse impulsionnelle finie qui varie d'un echantillon 
a 1' autre. II caracterise d' une maniere indirecte la 
realisation du canal multi-tra j ets durant la tranche 
temporelle de reception correspondante . 

Les echantillons en sortie du filtre blanchisseur 
correspondant a la sequence de reference sont utilises 
pour estimer le canal discret [1, 4, 5]. Cette 
estimation du canal discret permet d'egaliser les 
echantillons restants permettant eventuellement de 
detecter les deux sequences de donnees emises. 
L'egaliseur habituellement utilise est connu sous le 
nom d'egaliseur de VITERBI avec decisions souples [2]. 
Comme son nom I'indique, cet egaliseur fournit, grace a 
une application de I'algorithme de VITERBI avec sorties 
ponderees ("Soft-Output Viterbi Algorithm" ou SOVA en 
anglais) [2], des decisior^s souples sur tous les 
symboles de donnees emis. Ces sorties souples sont 
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eventuellement desent relacees et decodees pour detecter 
les informations envoyees . 

La figure 2 annexee illustre ces operations 
d ' emission/reception . L'emetteur E comprend une source 
de donnees 10, un codeur/entrelaceur 12 delivrant des 
symboles a^, un modulateur 14. Le canal multi-tra j ets 
est symbolise par le bloc 20. Le recepteur R comprend 
un demodulateur 30 (filtre adapte/f litre blanchisseur ) 
delivrant des echantillons Rk, un estimateur de canal 
discret 32, un egaliseur 34, un desentrelaceur/decodeur 
36 et finalement un destinataire 38. 

Le canal discret vu a la sortie du filtre 
blanchisseur peut varier de maniere significative d' une 
tranche temporelle a la suivante. Cette variation est 
due principalement au changement des conditions de 
propagation entre l'emetteur et le recepteur et a la 
stabilite en frequence au niveau du recepteur. 

Les conditions de propagation ont une influence 
directe sur le canal multi-tra j ets observe. Elles 
changent en raison d'une modification de 
1' environnement ou d ' un deplacement de l'emetteur et/ou 
du recepteur. Elles entrainent une variation temporelle 
du canal discret, aussi bien entre tranches temporelles 
successives qu'a I'interieur d'une meme tranche. 

La variation du canal discret entre deux tranches 
temporelles successives allouees au meme utilisateur 
est d'autant plus grande que I'intervalle temporel 
entre ces tranches est grand. Cette variation est 
accentuee, meme au niveau d'une seule tranche 
temporelle, en raison d'une augmentation de la 
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frequence porteuse ou de la vitesse de I'emetteur ou/et 
du recepteur. 

Dans la pratique, la variation du canal discret 
entre deux tranches temporelles est suf f isamment grande 
5 pour empecher toute estimation adaptative de celui-ci. 
L' estimation de canal au niveau d'une tranche 
temporelle doit alors reposer sur les seuls 
echantillons de la sequence de reference 
correspondante . 

10 L' utilisation de cette sequence de reference 

fournit une estimation non biaisee du canal discret au 
milieu de la tranche temporelle. A I'instar du systeme 
GSM, on fait generalement I'hypothese d'un canal 
invariant au cours d'une tranche temporelle regue. Dans 

15 ce cas precis, 1' estimation obtenue au milieu de la 
frame peut servir sans degradation a 1'' egalisation du 
reste de la trame. Mais, le deplacement a grande 
Vitesse de certains terminaux et la demande soutenue en 
services a hauts debits operant a des frequences 

20 radioelectriques de plus en plus grandes, rendent cette 
hypothese de moins en moins justifiee. Le canal discret 
peut en effet subir des variations signif icatives entre 
le debut et la fin d'une tranche temporelle donnee. 
L' ecart entre ce canal discret reel et son estimation 

25 se creuse de plus en plus lorsqu'on s ' eloigne de la 
sequence de reference. Ceci peut entrainer une 
degradation importante et irreversible de la qualite de 
la reception et/ou des performances du systeme AMRT . 

Une solution [13] a ^deja ete proposee pour 

30 resoudre ce probleme de variation rapide du canal dans 
le cas des canaux radioelectriques lies aux 
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radiocommunications par satellite. Cependant, cette 
solution ne traite que le cas d'un canal avec un seul 
trajet et suppose que les symboles, de reference ou de 
donnees, ont une energie transmise commune. Elle ne 
5 peut done pas etre appliquee aux canaux multi-tra j ets 
rencontres dans les communications radiomobiles 
terrestres . 



Par ailleurs, a I'instar de 1' extension EDGE du 





systeme 


GSM 


, des debits importants peuvent etre 


10 


obtenus , 


a bande f requentielle inchangee, par la mise 




en oeuvr 


e de 


modulations utilisant un plus grand nombre 




d' etats . 


Le 


systeme devient alors tres sensible a la 




qualite 


de 1' 


'estimation de canal, et 11 peut subir des 




degradat 


ions 


importantes meme a de tres faibles 


15 


vitesses 


de 


deplacement, si la puissance et/ou la 




longueur 


de 


la sequence de reference ne sont pas 



augmentees . 

Le but de la presente invention justement est de 
remedier a ces inconvenients . 

20 

Expose de 1 ' invention 

La presente invention a pour but d'accroltre les 
performances des systemes AMRT en reduisant le rapport 
entre la puissance regue et 1' interference d' acces 

25 multiple pour une qualite de reception donnee ou, en 
d'autres termes, d'ameliorer la qualite de reception a 
puissance de reception constante. Cette amelioration de 
la qualite permet, entre autres, d' augmenter la 
capacite et la couverture d^s systemes AMRT, tels que 

30 le systeme GSM et surtout son extension EDGE pour 
offrir des services a hauts debits. Cette amelioration 
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est obtenue notamment grace a une optimisation de 
1' estimation de canal au niveau du recepteur. Cette 
optimisation permet de contrecarrer toute degradation 
en performances engendree par une variation rapide du 
5 canal entre le debut et la fin d'une meme tranche 
temporelle. Elle permet egalement de reduire 
signif icativement les degradations engendrees par une 
augmentation du nombre d'etats de la modulation et ceci 
sans recours a une augmentation de la longueur ou de la 

10 puissance de la sequence de reference. 

Un autre but de 1' invention est de reduire le prix 
de revient des terminaux en les rendant, a qualite de 
reception inchangee, moins sensibles a 1 ' instabilite en 
frequence d' un oscillateur local bon marche permettant, 

15 entre autres, la transposition du signal regu en bande 
de base. 

L' invention permet encore de reduire, a qualite de 
reception egale, le nombre et/ou la puissance des 
symboles de reference transmis durant chaque tranche 

20 temporelle. Ce but est atteint par la prise en compte, 
dans 1' estimation du canal discret, des symboles de 
donnees qui sont en general plus nombreux que les 
symboles de reference. Ce but est egalement atteint par 
la prise en compte de maniere optimale des echantillons 

25 correspondant aux symboles de reference et/ou de 
donnees d'un nombre arbitraire de tranches temporelles 
allouees a 1 ' utilisateur . Cette prise en compte est 
effective meme en cas de variations importantes du 
canal entre deux tranches successives appartenant au 

30 meme utilisateur. Ce but est enfin atteint en prenant 
en compte de maniere optimale les echantillons 
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provenant du canal de signalisation commun. et 
d' eventuels utilisateurs sur la liaison descendante. 

Selon 1' invention, les symboles de reference sont 
repartis en groupes sur toute la duree de la tranche 
temporelle allouee, ce qui garantit une meilleure 
poursuite de la variation du canal et done une 
robustesse accrue vis a vis des grandes vitesses de 
deplacement et de 1 ' instabilite en frequence de 
1 ' oscillateur local. Les sequences de reference peuvent 
etre choisies independamment les unes des autres, mais 
elles doivent toujours presenter de bonnes proprietes 
d' autocorrelation . Elles peuvent etre derivees de 
sequences de type CAZAC de longueur reduite . 

Avec 1' invention, I'energie transmise par symbole 
peut varier d' un symbole a 1' autre sans que cela pose 
des difficultes pour 1' estimation du canal, laquelle 
s'effectuera toujours de maniere optimale. 

Selon 1' invention, le traitement d'une tranche 
temporelle se fait par bloc a chaque fois que les 
echantillons correspondant a cette tranche et 
eventuellement a ceux d' autres tranches (provenant du 
meme utilisateur, d' autres utilisateurs ou du canal de 
signalisation commun) sont disponibles pour effectuer 
1' estimation de canal. 

A I'instar d' un recepteur AMRT classique (tel que 
celui du systeme GSM) on commence toujours par une 
estimation grossiere du canal discret en utilisant 
uniquement des echantillons dependant des seuls 
symboles de reference de ^la tranche temporelle a 
traiter. Cette estimation grossiere permet de 
caracteriser, au niveau de chaque symbole de reference 
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ou de donnee du bloc a traiter, symbols par symbole, 
1' evolution de la phase et de 1' amplitude de tous les 
coefficients (egalement appeles "trajets" par abus de 
langage) du canal discret retenus pour la realisation 
5 de 1'' egalisation . 

L'une des caracteristiques de 1' invention est 
1 ' utilisation d'un treillis pour representer le canal 
multitra j ets . Cette caracteristique reduit la 
complexite du recepteur. En effet, si 1 ' on prend comme 

10 exemple des sequences de symboles de donnees comprenant 
n symboles, sans codage, avec une modulation de phase a 
deux etats (MDP2), chaque bit emis peut prendre l'une 
ou 1' autre de deux valeurs. Comme ces valeurs sont 
equiprobables, cela signifie qu'il existe 2" sequences 

15 de symboles de donnees possibles. La complexite d'un 
calcul traditionnel utilisant les probabilites 
partielles pour toutes les sequences possibles croit 
done exponentiellement (2") avec le nombre de symboles 
par sequence. 

20 L ' utilisation, conformement a 1' invention, d'un 

treillis pour representer le canal reduit 
considerablement la complexite des calculs, car ceux-ci 
vont se subdiviser en sections, chaque section donnant 
lieu a un certain type de calcul, toujours le meme pour 

25 chaque section. La complexite augmente done comme le 
nombre de sections, lequel croit comme le nombre n. La 
complexite n'est done plus exponentielle mais lineaire. 

Les etiquettes des branches des sections du 
treillis changent au gre de ^ la variation temporelle du 

30 canal discret au niveau de chacune des tranches 
temporelles constitutives du bloc. Par ailleurs, un 
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algorithme different de 1' algorithms de VITERBI .avec 
sorties souples (SOVA) [2] est utilise pour une 
caracterisation plus precise des probabilites des 
branches du treillis conditionnellement aux 
5 echantillons regus. Dans le cas d'une mise en oeuvre 
optimale de I'estimateur de canal, un algorithme dit de 
BAHL ou BCJR (abreviation de "Bahl Cocke Jelinek 
Raviv") [6] peut etre utilise pour fournir ces 
probabilites conditionnelles . D'autres algorithraes [10] 

10 plus simples, provenant generalement d'une 
linearisation de cet algorithme, peuvent etre utilises, 
avec cependant une legere degradation en performance, 
lis permettent de reduire considerablement la 
complexite de I'estimateur de canal. 

15 Dans le cas d'un recepteur classique, les seules 

sorties ponderees fournies par 1' algorithme SOVA sont 
directement exploitees apres desentrelacement pour 
decoder les informations transmises. Dans le cas du 
recepteur de 1' invention, les probabilites 

20 conditionnelles caracterisant indirectement les 
symboles de donnees sont exploitees, en plus des 
symboles de reference, pour fournir une estimation de 
meilleure qualite du canal discret. 

L' estimation amelioree du canal discret obtenue a 

25 la fin d'une iteration donnee peut etre utilisee pour 
augmenter la qualite et la fiabilite des probabilites 
conditionnelles des branches du treillis caracterisant 
le canal discret. Ces probabilites conditionnelles 
ameliorees peuvent, a leur tour, etre utilisees, 

30 conj ointement avec les symboles de reference, pour 
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apporter une amelioration supplementaire a 1' estimation 
du canal discret. 

En theorie, une estimation iterative ne conduit a 
une estimation optimale du canal discret qu'apres un 
5 nombre infini d ' iterations . En pratique cependant, 
quelques iterations suffisent sans grande degradation 
des performances par rapport a une estimation optimale. 

Une fois 1' estimation du canal discret accomplie, 
le treillis correspondant peut etre utilise soit par 
10 I'algorithme optimal de BAHL (ou tout autre 
simplification de celui-ci) , soit par I'algorithme SOVA 
pour fournir des sorties ponderees ameliorees au 
decodeur . 

Le caractere optimal du recepteur selon 

15 1' invention est lie a la qualite de 1' estimation du 
canal discret vu durant une tranche temporelle donnee. 
Ce caractere optimal repose sur 1' utilisation d'un 
algorithme iteratif dit SAGE (abreviation de "Space- 
Alternating Generalized Expectation-Maximization 

20 Algorithm") [7] pour trouver la realisation de canal la 
plus vraisemblable condit ionnellement aux echantillons 
regus du bloc a traiter. L' algorithme SAGE est une 
extension de I'algorithme EM (contraction de 
"Expectation-Maximization") [8, 9] qui permet de se 

25 debarrasser du probleme de couplage entre coefficients 
du canal discret lors de leur estimation. 

L' estimation du canal discret repose egalement sur 
la decomposition de chaque trajet du canal discret a 
estimer selon un algorit^hme d' expansion dit de 

30 KARHUNEN-LOEVE [3]. Cette decomposition fournit, d' une 
part, une caracterisation souple des variations 
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temporelles de chacun des trajets discrets et, d'autre 
part, s'integre naturellement dans I'algorithme SAGE. 
Elle permet egalement de traiter le cas ou I'energie 
emise par symbole varie d' un symbole a 1' autre. 
5 L' estimation du canal discret repose enfin sur 

1' utilisation de I'algorithme de BAHL, qui fournit a 
chaque iteration de I'algorithme SAGE les probabilites 
des branches du treillis caracterisant le canal 
discret, conditionnellement aux echantillons regus et a 
10 1' estimation de canal fournie a 1' iteration precedente . 



De fagon precise, 1' invention a done pour objet un 
precede de communications radiomobiles du type a acces 
multiple a repartition en temps, dans lequel on alloue 
15 a plusieurs utilisateurs des tranches temporelles pour 
1' emission et la reception de signaux radioelectriques, 
et dans lequel : 

- a 1' emission, on constitue, dans chaque tranche 
temporelle allouee, des symboles de donnees et 

20 des symboles de reference, on emet un signal 

radioelectrique contenant ces symboles de 
donnees et ces symboles de reference, ce signal 
empruntant un canal radiomobile mult i-tra j ets , 

- a la reception, on effectue un filtrage adapte 
25 du signal regu, on echantillonne le signal 

filtre au rythme des symboles, on effectue une 
estimation du canal radiomobile, on traite les 
echantillons en tenant compte de cette 
estimation du canal et 1 ' on restitue les 
30 symboles de donnees propres a la tranche 

allouee. 



SP 16348. C/RS 



13 

precede etant caracterise en ce que : 

- a 1' emission, pour chaque tranche temporelle 
allouee, on constitue plusieurs sequences de 
reference possedant des proprietes de correlation 
temporelle appropriees et 1 ' on repartit ces 
sequences dans chaque tranche temporelle, 

- a la reception : 

i) pour estimer le canal radiomobile, on 
represente ce canal par un treillis 
comprenant des branches symbolisant les 
transitions entre deux etats successifs du 
canal et caracterisant les sequences 
possibles a la sortie du canal, 
ii) on calcule les probabilites de ces branches 
en mettant en ceuvre un algorithme iteratif 
dans lequel on calcule, par un nombre fini 
(D+1) d' iterations et selon le critere de 
probabilite maximum a posteriori, une 

estimation d'une representation [G^^'^^jf^Q des 
L trajets de ce canal, 1 ' iteration de rang 
(d+1) (3d+l) permettant d'obtenir une 

reestimation {c^^'^^'^j^^o a partir de 
1' estimation {g^^'^^}J'J obtenue par 1' iteration 
precedente (3d), et, a chaque iteration 
(3d+l) , on calcule en L etapes les vecteurs 
jQf(d+i)j^ ^=0, 1, L-1 a partir des 

vecteurs {g^^^^^^^ et {G^^^'^^f, calcules a 
I'etape precedente, 1' estimation des 
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L trajets du canal etant finalement obtenue 
a partir de 1 ' est imation de sa 

representation {g^^^^I^Jo • 

Breve description des dessins 

- la figure 1, deja decrite, montre la repartition 
des tranches teraporelles allouees a plusieurs 
utilisateurs ; 

- la figure 2, deja decrite, est un schema 
illustrant 1' emission et la reception dans la 
technique AMRT ; 

- la figure 3 montre les variations du taux 
d'erreur binaire (TEB) en fonction de 1 ' etalement 
Doppler normalise (BdTs) ; 

- la figure 4 montre les variations du taux 
d'erreur binaire (TEB) en fonction du rapport 
Eg/No pour des modulations MDP4 et MDP8, dans une 
premiere configuration et pour un canal discret a 
deux trajets ; 

- la figure 5 montre les variations du taux 
d'erreur binaire (TEB) en fonction du rapport 
Eg/No pour des modulations MDP4 et MDP8, dans une 
deuxieme configuration avec un canal discret a 
deux trajets ; 

- la figure 6 montre les variations du taux 
d'erreur binaire (TEB) en fonction du rapport 
Eg/No pour une modulation MDP4 utilisant la 
deuxieme configuration et un canal discret a 
quatre trajets ; 
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- la figure 7 montre les variations du "taux 
d'erreur binaire (TEB) en fonction du rapport 
Eg/N pour une modulation MDP4 utilisant un canal 
discret avec un, deux et quatre trajets ; 

5 - la figure 8 est une representation du canal 

discret vu a la sortie du filtre blanchisseur, 
sous forme de registre a decalage ; 

- la figure 9 illustre 1' estimation des 
coefficients du filtre du canal discret par 

10 correlation des echantillons regus avec une 

sequence CAZAC et ses decalees circulaires ; 

- la figure 10 montre la repartition des valeurs 
propres les plus importantes de la matrice de 
covariance pour trois etalements Doppler 

15 normalises ; 

- la figure 11 illustre une representation en 
treillis du canal discret normalise vu a la 
sortie du filtre blanchisseur pour un canal a 
deux trajets et une modulation a deux etats et 

20 une longueur de tranche temporelle egale a 

quatre ; 

- la figure 12 montre la chronologie des 
traitements effectues a la reception pour une 
tranche temporelle regue ; 

25 - la figure 13 montre la chronologie des 

traitements effectues dans une iteration 
elementaire de 1 ' algorithme d' estimation de 
canal. y 
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Description de modes particuliers de mise en ceuvre • 

Les performances du precede selon 1' invention ont 
ete comparees a celles d'une variante du precede 
classique utilise dans le cadre du systeme GSM [l].Le 
5 recepteur classique du systeme GSM estime d' abord le 
canal discret au milieu de la tranche temporelle a 
traiter, par correlation des echantillons regus 
dependant uniquement de la sequence de reference avec 
la sequence CAZAC correspondante . II egalise ensuite 

10 les autres echantillons regus au moyen de I'algorithme 
SOVA afin de decoder les symboles d' information emis . 
Afin d' evaluer equitablement 1' amelioration intrinseque 
de performances apportee par 1' invention par rapport au 
recepteur classique du GSM, I'algorithme sous-optimal 

15 SOVA a ete remplace par I'algorithme optimal de BAHL 
utilise dans 1' invention. 

La comparaison est obtenue a travers 1' evolution 
du taux d'erreur binaire (TEB) brut (sans prise en 
compte d'un eventuel codage correcteur ou detecteur 

20 d'erreurs), pour des symboles de reference aussi bien 
groupes au milieu de la tranche temporelle a traiter (a 
la maniere du systeme GSM) que repartis en petits 
groupes sur toute la duree de cette tranche 
(conf ormement a 1 ' invention) , afin de garantir une 

25 meilleure poursuite de la variation temporelle du 
canal. 

La figure 3 represente ainsi 1' evolution du TEB en 
fonction de I'etalement Doppler normalise B^T^, produit 
de I'etalement Doppler 5^ et de la duree d'un symbole 
30 , qui caracterise I'intensite de variation du canal 

due au deplacement de I'emetteur et/ou du recepteur. 
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L'etalement Doppler est a la fois proport ionnel "a la 
Vitesse du terminal mobile et a la frequence de la 
porteuse. L ' instabilite en frequence de I'oscillateur 
local, qui n'est pas prise en compte ici, conduit a des 
5 degradations supplementaires en performances. A 
I'instar du systeme GSM, le comportement du TEB a ete 
evalue pour des tranches temporelles composees de N=148 
symboles chacune. La modulation retenue est la 
modulation MDP4 (modulation a deplacement de phase a 4 
10 etats) . Le canal discret considere est compose de deux 
trajets (L=2) de meme puissance moyenne. Le rapport 
Eg/M entre I'energie moyenne regue par symbole et le 
niveau de bruit est de 10 dB. Trois configurations ont 
ete considerees : 
15 ■ Une sequence de reference centrale de 26 symboles 

bipolaires prise au hasard parmi les 8 sequences 
utilisees dans le systeme GSM. Ces sequences sont 
des extensions cycliques des 8 sequences CAZAC de 
longueur 16 [1, 11] . 
20 ■ Une sequence de reference centrale, {-1, +1, +1, 

+ 1, -1}, de 5 symboles bipolaires obtenue par 
extension cyclique (a un signe et un decalage 
pres) de 1' unique sequence CAZAC bipolaire de 
longueur 4 a savoir {-1, +1, +1, +1}. 
25 ■ Cinq groupes de sequences de reference, repartis 

sur toute la tranche temporelle a traiter, 
identiques a la sequence de reference centrale de 
5 symboles bipolaires precedents (25 symboles de 
reference au total). ' 

30 
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Sur la figure 3, les courbes representees 
correspondent aux cas suivants : 

41 : recepteur classique, premiere configuration, 

42 : recepteur classique, deuxieme configuration, 

43 : recepteur de 1' invention, deuxieme 
configuration, 

44 : recepteur classique, troisieme 
configuration, 

45 : recepteur de 1' invention, troisieme 
configuration . 

Ces resultats montrent que 1' utilisation de 
sequences groupees au centre de chaque tranche 
temporelle (premiere et deuxieme configurations) ne 
permet d' aller au-dela d' un etalement Doppler normalise 
BdTs=1/2500 sans degradation significative des 
performances. A I'inverse, 1 ' utilisation de-groupes de 
petites sequences reparties (troisieme configuration) 
permet d'atteindre, sans degradation perceptible, des 
etalements Doppler normalises proches de BdTs=1/100. 
Cette amelioration significative de la robustesse du 
recepteur vis-a-vis de I'effet Doppler, et 
eventuellement de 1 ' instabilite en frequence de 
1' oscillateur local, est due a la repartition des 
sequences de reference qui permet de suivre de maniere 
plus precise toute evolution temporelle du canal. 

Le tableau 1 caracterise de maniere chiffree 
1' amelioration en performances apportee par 1' invention 
vis-a-vis de I'algorithme cla^ssique ameliore du systeme 
GSM. Dans ce tableau R.C. signifie "recepteur 
classique" et R.I. recepteur selon 1' invention. Le 
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rapport Eg/NO entre 1 ' energie moyenne regue par symbols 
et le niveau de bruit est de 10 dB. Dans le cas d'une 
sequence de reference centrale longue (premiere 
configuration), 1 ' amelioration apportee par 1' invention 
5 est negligeable. 



Etalement 
DOPPLER 
normalise 


Premiere 
configuration 


Deuxieme 
configuration 


Troisieme 
configuration 


R.C. 


R.I. 


R.C. 


R.I. 


R.C. 


R.I. 


1/2500 


0,00 
99 


0,00855 


0,01584 


0,0089 


0,0105 


0,0077 


1/100 


X 


X 


x 


X 


0,0143 


0,0083 



Tableau 1 

10 

La sequence de reference, dans ce cas, est 
suffisamment longue pour garantir, sans I'aide des 
sequences de donnees inconnues, una bonne qualite 
d' estimation de canal. Cependant, dans le cas d'une 

15 seule sequence de reference courte (deuxieme 
configuration), le recepteur selon 1' invention permet 
d'obtenir non seulement de meilleures performances que 
le recepteur classique, mais aussi des performances 
comparables a ceux obtenus avec la sequence de 

20 reference longue de 26 symboles de reference. Par 
consequent, le recepteur selon 1' invention est rendu, 
d'une certaine maniere, insensible a une reduction meme 
importante de la taille de la sequence de reference. Ce 
resultat s'explique par la prise en compte, dans 

25 1' estimation de canal, des 'symboles des sequences de 
donnees qui viennent cooperer avec les symboles de 
reference. Dans le cas de cinq groupes de sequences de 
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reference courtes (troisieme configuration) , • le 
recepteur selon 1' invention permet de maintenir le taux 
d'erreur binaire presque inchange en faisant varier 
I'etalement Doppler normalise entre 1/2500 et 1/100. 

5 

Les figures 4, 5 et 6 representent le taux 
d'erreur binaire en fonction du rapport E^fN^ entre 
I'energie raoyenne regue par symbole et le niveau du 
bruit pour un etalement Doppler normalise BdTs=1/2500 

10 peu contraignant . 

Pareillement au systeme GSM, le comportement du 
TEB a ete evalue pour des tranches temporelles 
composees de N==148 symboles chacune. Toutefois, les 
modulations retenues sont les modulations MDP4 et MDP8 . 

15 Le canal discret considere est compose de un (L=l) ou 
de deux (L=2) ou de quatre (L=4) trajets de meme 
puissance moyenne. Deux configurations ont ete 
considerees : 

■ Une sequence de reference centrale de 26 symboles 
20 bipolaires prise parmi les 8 sequences utilisees 

dans le systeme GSM. Ces sequences sont des 
extensions cycliques des 8 sequences CAZAC de 
longueur 16 [1, 11] . 

■ Une sequence de reference centrale ({-1, +1, +1, 
25 +1, -1} pour deux trajets et {-1, +1, +1, +1, -1, 

+1, +1} pour quatre trajets de 4+L-l symboles 
bipolaires, obtenue par extension cyclique de 
1' unique (a un signe et un decalage pres) 
sequence CAZAC bipolaire de longueur 4 a savoir 
30 {-1, +1, +1, +1}. 
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Pour la figure 4, les modulations retenues sont- les 
modulations MDP4 et MDP8 utilisant la premiere 
configuration avec deux trajets (L=2). Les courbes se 
rapportent aux cas suivants : 
5 - 51 : MDP8, recepteur classique, 

- 52 : MDP8, recepteur selon 1' invention, 

- 53 : MDP8, recepteur avec canal discret 

parfaitement connu, 

- 54 : MDP4, recepteur classique, 

10 - 55 : MDP4, recepteur selon 1' invention, 

- 56 : MDP4, recepteur avec canal discret 

parfaitement connu. 

Sur la figure 5, les modulations retenues sont les 
15 modulations MDP4 et MDP8 utilisant la deuxieme 
configuration et un canal discret avec deux trajets 
(L=2) . Les courbes se rapportent aux cas suivants : 

- 61 : MDP8, recepteur classique, 

- 62 : MDPB, recepteur de 1' invention, 

20 - 63 : MDPB, recepteur avec canal discret 

parfaitement connu, 

- 64 : MDP4, recepteur classique, 

- 65 : MDP4, recepteur de 1' invention, 

- 66 : MDP4, recepteur avec canal discret 
25 parfaitement connu. 



Pour la figure 6, la modulation retenue est la 

modulation MDP4 avec la deuxieme configuration et un 

canal discret avec quatre trajets (L=4). Les courbes se 
30 rapportent aux cas suivants : 
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- 71 : recepteur 

- 72 : recepteur 

- 73 : recepteur 

connu . 

5 

Les courbes des figures 4 et 5 montrent que, a la 
difference d'un recepteur selon 1' invention, les 
performances du recepteur classique du systeme GSM sont 
tres sensibles a la longueur de la sequence de 

10 reference utilisee et se degradent substantiellement 
quand on passe de la premiere configuration a la 
deuxieme. Les courbes des figures 5 et 6 montrent, en 
outre, que les degradations en performances du 
recepteur classique augmentent signif icativement avec 

15 le nombre de trajets alors que, a 1' inverse, ces 
degradations restent presque inchangees et faibles pour 
un recepteur selon 1' invention. Le recepteur selon 
1' invention presents done une excellente robustesse vis 
a vis de toute reduction de la longueur de la sequence 

20 de reference et de toute augmentation du nombre de 
trajets du canal. 

Le tableau 2 caracterise de maniere chiffree, dans 
le cas de la modulation MDP4, 1' amelioration en 
25 performances apportee par le recepteur selon 
1' invention vis-a-vis de I'algorithme ameliore du 
systeme GSM. Le rapport BqTs est egal a 1/2500, le taux 
d'erreur binaire est egal a 10"^ et N=148. 



classique, 
de 1' invention, 
avec canal discret parfaitement 
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1 Gain en dB 


Nombre de 


Sequence de 


Sequence de 


Recepteur selon 1' invention avec 


trajets {L) 


reference 


reference 


sequence courte, recepteur 




longue 


courte 


classique avec sequence longue 


1 


0,03 


1,03 


0,03 


2 


0,48 


1,83 


0,48 


4 


1,42 


3,04 


1,43 



Tableau 2 



5 Ce tableau montre que, pour une longueur donnee 

(premiere ou deuxieme configuration) de la sequence de 
reference, un recepteur selon 1' invention a toujours de 
meilleures performances qu'un recepteur classique 
ameliore. Ces resultats montrent aussi que 

10 1' amelioration apportee par 1' invention est d'autant 
plus grande que la longueur de la sequence de reference 
est reduite ou que le nombre de trajets du canal est 
augmente. Ces resultats montrent enfin que, meme avec 
une sequence de reference courte (deuxieme 

15 configuration), un recepteur selon 1' invention a 
toujours de meilleures performances qu'un recepteur 
classique avec une sequence de reference longue 
(premiere configuration) . 

20 La figure 7 represente le taux d'erreur binaire en 

fonction du rapport / Nq pour un etalement Doppler 
normalise BdTs=1/100 contraignant et des tranches 
temporelles composees de N=100 symboles chacune . La 
modulation retenue est toujours la modulation MDP4 . Le 

25 canal discret considere est compose de un (L=l) ou de 
deux (L=2) ou de quatre (L=4) trajets de meme puissance 
moyenne . Une seule configuration avec trois sequences 
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de reference reparties est consideree. Chaque sequence 
est composee de 4+L-l symboles bipolaires et est 
obtenue par extension cyclique (a un signe et un 
decalage pres) de 1' unique sequence CAZAC bipolaire de 
longueur 4 a savoir {-1, +1, +1, +1}. 

La valeur de 1 ' etalement Doppler consideree est 
suf f isamment contraignante pour ne pas permettre a un 
recepteur classique de suivre les variations de chacun 
des trajets du canal discret qui peuvent subir des 
rotations de phase de ±7r entre le debut et la fin de 
chaque tranche temporelle. 

Les courbes de la figure 7 se rapportent aux cas 
suivants : 

- 81 : L=l, recepteur avec canal discret 

parf aitement connu, 

- 82 : L=l, recepteur de 1' invention, 

- 83 : L=2, recepteur avec canal discret 

parf aitement connu, 

- 84 : L=2, recepteur de 1' invention 

- 85 : L=4, recepteur avec canal discret 

parf aitement connu, 

- 86 : L=4, recepteur de 1' invention. 

Les courbes de la figure 7 montrent que le 
recepteur selon 1' invention conduit a une degradation 
tres faible des performances vis-a-vis d' un recepteur 
ideal connaissant parf aitement le canal discret. Cette 
degradation est de I'ordre de 0,4 dB pour un seul 
trajet (L=l), 0,5 dB pour deux trajets (L=2) et 1 dB 
pour quatre trajets (L=4). 
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Les performances du precede de 1' invention ayant 
ete exposees, la phase de reception et notamment le 
traitement du signal regu, vont etre decrits plus en 
details . 

Un recepteur selon 1' invention determine, 
conditionnellement aux echantillons regus du bloc a 
traiter, la realisation du canal discret la plus 
vraisemblable durant une tranche temporelle donnee. 
Cette determination se fait selon un critere dit de 
probabilite a posteriori maximale (MAP en abrege pour 
"Maximum A Posteriori probability" en anglais) . Cette 
estimation du canal necessite la connaissance des 
probabilites a priori des symboles des sequences de 
donnees transmises. 

Par souci de simplification, les developpements 
mathematiques qui suivent et les figures 
correspondantes se rapportent toujours au cas ou chaque 
bloc a traiter pour 1' estimation de canal est compose 
uniquement des echantillons regus se rapportant a une 
seule tranche temporelle. Le cas ou les symboles (de 
reference ou non) d' autres tranches teraporelles (du 
meme utilisateur ou d' autres ut ilisateurs , du canal 
comm.un de signalisat ion, . . . ) sont utilises se traite 
de la meme fagon. 

On note I'operateur de transposition. Comme 

indique en figure 2, toute tranche temporelle d' un 
utilisateur donne est composee d'une suite de N 
symboles aQ,a,,. . .,a^_[ que 1 ' on peut ecrire sous forme de 
vecteur colonne 
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les symboles etant eventuellement modules en phase 
(MDP2, MDP4, MDP8, ...)- L'energie transmise E,^ 
associee a chaque symbole a,^ est independante des 
valeurs prises par ce symbole mais peut varier d' un 
5 symbole de la tranche temporelle a 1' autre. Cette 
variation est generalement la consequence d' une 
allocation par le systeme AMRT d' energies differentes 
aux sequences de donnees et aux sequences de reference. 
Comme indique egalement dans la meme figure, les 
10 echantillons en sortie du filtre blanchisseur 
correspondant a cette tranche temporelle forment un 
vecteur colonne 

de A^+Z-l composantes, ou L est le nombre de trajets 

15 du canal discret. 

La figure 8a rassemble les moyens compris entre 
1' entree du modulateur de I'emetteur et la sortie du 
filtre blanchisseur du recepteur. Cette figure fournit 
done une representation, du canal radioelectrique reel 

20 et du traitement effectue au niveau du recepteur 
(demodulation, filtrage adapte et filtrage 

blanchissant ) . Sur cette figure, on trouve le 
modulateur 91, un bloc 92 symbolisant le canal multi- 
trajets, le demodulateur 93 avec son filtre adapte 94 

25 et son filtre blanchisseur 95. 

La figure Bb montre une representation de ces 
moyens sous forme de registre a decalage avec (L-l) 
cellules Yi, Y2, Yl-i- Ce registre regoit le symbole 

courant an et les cellules contiennent les symboles 

30 retardes d'une periode symbole, a savoir an-i, stn-z, ••■r 
an-(L-i)- Le circuit represente comprend encore (L) 
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multiplieurs Zo, Zi, . . . , Zl-i recevant respectivement 
(L) coefficients cj, cj,, ... cj~' , un addit ionneur S 
relie aux multiplieurs, un dernier additionneur 
recevant un signal Wn representant un bruit. Get ultime 
additionneur delivre un signal rn- 

Chaque echantillon regu peut se mettre sous la 

forme 

R,=Y^ia,_i+W, , k = 0X...,N + L-2 , 

1=0 

ou cl est le /erne coefficient du filtre de canal 
discret au A: erne instant d' echantillonnage et est un 
bruit additif blanc gaussien. Pour simplifier la 
presentation et avoir des equations concises, les 
coefficients du canal discret, verifiant k<l ou 

k>N + l, interviennent dans 1' expression de R,, alors 
qu' ils sont nuls. Pour la meme raison, les symboles 
emis, verifiant A; < 0 ou k>N-l, apparaissent dans 
cette expression de i?^. meme s'ils sont nuls. 

Les coefficients du canal discret sont presque 
decorreles entre trajets discrets differents 
(independance vis-a-vis de la variable /). Abstraction 
faite des effets de la modulation et de la 
demodulation, cette propriete est due a 1' independance 
des trajets multiples du canal radioelectrique reel. 
Toutefois, les coefficients d'un meme trajet discret / 
sont correles entre eux (dependance vis-a-vis de la 
variable k) . Leur fonction d' autocorrelation discrete, 
definie par 
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depend uniquement du spectre de puissance Doppler du 
trajet correspondant . 

Les spectres de puissance Doppler generalement 
rencontres en communications radiomobiles sent soit de 
type classique^ (en environnements exterieurs) , soit de 
type plat, (en environnements interieurs). La fonction 
d' autocorrelation du I erne trajet, de puissance moyenne 
4' , est donnee par 

dans le cas d' un spectre de puissance Doppler classique 
et par 

dans le cas d' un spectre de puissance Doppler plat, ou 
Jo(-) est la fonction de Bessel de premiere espece 
d'ordre 0 [3, 12]. La puissance moyenne 4'1> varie d' un 
trajet a 1' autre et caracterise le profil de 
I'intensite des trajets ("multipath intensity profile" 
en anglais [3]) du canal discret. 

Le recepteur de 1' invention doit determiner selon 
le critere du maximum a posteriori (MAP) , la 
realisation la plus vraisemblable des vecteurs 

c' =(c;,c;,i,c;,2,...,c;,,_5)' , / = o,i,...,z-i , 

caracterisant les L trajets discrets retenus par le 
recepteur pour 1 ' egalisation . A cette fin, il a besoin 
d'une representation adequate du canal discret lui 
permettant de s'affranchir de la correlation temporelle 
des coefficients d' un meme trajet. II a egalement 
besoin d'une representatio/n qui lui permette de 
s'affranchir de la variation eventuelle de I'energie 
transmise d' un symbole de la tranche temporelle a 
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1' autre. Cette representation est obtenue tout d'^abord 
par 1' introduction du vecteur de symboles transmis 
normalises 

5 avec des composantes A^=aJ.[EldQ module constant egal 
a un. Elle est egalement obtenue en reecrivant les 
echantillons regus durant une tranche temporelle sous 
la forme : 

L-\ 

10 ou C[ est la keme composante du vecteur normalise 

associe et caracterisant directement le /erne trajet du 
canal discret. Elle est enfin acquise grace a une 
decomposition, dite de KARHUNEN-LOEVE [3], de chaque 

15 vecteur normalise C' dans une base orthonormale {B*}^^^ 
tenant compte a la fois de la variation de I'energie 
des symboles emis et des caracterist iques du spectre de 
puissance Doppler du /erne trajet [3]. Cette 
decomposition permet en effet de decrire le vecteur 

20 normalise C' (avec composantes temporelles correlees) 
du /erne trajet discret par un vecteur colonne 
equivalent, G' = (G^,G;,...,G;,_,)' , dont les composantes 
sont bien decorrelees. Le vecteur C' peut etre 
reconstruit a partir de sa representation adequate G 

25 par la transformation 
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Inversement, les composantes de cette 

representation adequate peuvent etre retrouvees a 

partir du meme vecteur C' par la transformation inverse 

=(B'*)''C' , ;t = 0,1,...,A^-1, 
ou (•)' est I'operateur de conjugaison. 

Conf ormement aux proprietes de la decomposition de 
KARHUNEN-LOEVE [3], les vecteurs de la base 
orthonormale associee au /erne trajet ne sont autres que 
les vecteurs propres de la raatrice de covariance H du 
vecteur normalise C' associe a ce trajet. De plus, les 
valeurs propres {-^*}^_o associes ne sont autres que les 
variances des composantes Go,G] ,...,G^_, de la 
representation adequate G' du trajet correspondant . Les 
vecteurs {b'*}^J"q associes a chaque trajet, ainsi que les 
valeurs propres correspondantes, peuvent done etre 
obtenus a partir de 1' equations matricielle 
H'B''=A"B*, k = OX...,N-\. 

La (m,«)eme entree de la matrice de covariance 
du vecteur normalise associe au I erne trajet est donnee 
directement par 

Hi„ = e[cUc'S] = ^E„EJ„,.„ , 

ou (p'^_„ denote la fonction d' autocorrelation discrete de 
ce trajet definie plus haut . 

Le recepteur selon 1' invention doit normalement 
disposer de L bases orthonormales associees a chacun 
des L trajets du canal discret afin de reconstruire ou 
decomposer toute estimation intermediaire ou finale du 
canal discret. Dans la pratique, les formes des 
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spectres de puissance Doppler des L trajets du canal 



associees sont done identiques, a un facteur 
multiplicatif pres, et les bases orthonormales 
correspondantes le sont aussi. Une seule base 
orthonormale est alors generalement suffisante au 
niveau du recepteur pour effectuer convenablement 
l' estimation du canal. 

Pour generer la base orthonormale du /erne trajet, 
le recepteur doit disposer d' une estimation aussi 
precise que possible de la matrice de covariance H' de 
ce trajet. A cette fin, il doit disposer d'une 
connaissance aussi precise que possible des 
caracteristiques du spectre de puissance Doppler du 
canal radioelectrique . Cette connaissance est mise a 
jour regulierement et comporte obligatoirement une 
estimation de I'etalement Doppler Bd ou d'une borne 
superieure de celui-ci. Elle peut comporter 
eventuellement une estimation de la puissance moyenne 
de chacun des trajets du canal discret et de la forme 
du spectre de puissance Doppler. Condit ionnellement a 
cette connaissance, le recepteur peut, soit calculer en 
temps reel les bases orthonormales adequates, soit 
rechercher dans une banque de bases orthonormales pre- 
calculees celles qui concordent le mieux avec les 
caracteristiques du canal reel. En general, le 
recepteur n' a ni connaissance precise de la forme du 
spectre de puissance Doppler, ni connaissance parfaite 
et instantanee de I'etalement Doppler. Dans ce cas 
typique, il peut eventuellement faire I'hypothese d'un 
spectre de puissance Doppler plat avec, comme etalement 



sont identiques. Les matrices de covariance 
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Doppler, une borne superieure de I'etalement Doppler 
reel. Par ailleurs, le recepteur peut ne pas disposer 
d'une estimation precise de la puissance moyenne de 
chacun des trajets du canal discret pris 
individuellement . Dans ce cas aussi, 11 peut aussi 
faire I'hypothese d' un canal discret avec L trajets de 
puissances moyennes egales. 

L' estimation des vecteurs c' caracterisant les L 
trajets du canal discret se fait iterativement selon le 
critere du MAP a travers les vecteurs normalises C ou 
leur representation adequate . Cette estimation peut 
prendre en compte de maniere optimale aussi bien les 
caracteristiques des symboles de donnees et de 
reference emis ( statistiques et energies) que les 
caracteristiques du canal multi-tra j ets . 

L' estimation au MAP {g'}^^J"J d'une realisation de la 
representation adequate {g'}^^J du canal discret est la 
valeur 

{g ' } = arg maxpi{G'} "jj iR) 
qui maximise la densite de probabilite conditionnelle a 
posteriori p({g'}|'_'J|R) . Conformement a 1' invention, et 
en utilisant 1' extension dite SAGE [7] de I'algorithme 
EM [8, 9], il est possible d'obtenir iterativement une 
solution aussi proche qu' il le faut de la solution 
exacte . 

En general, la densite de probabilite 
conditionnelle a posteriori p({g'}';Jr) a maximiser 
possede plusieurs maxima globaux conduisant souvent a 
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une ambiguite de 1 ' estimation de canal selon le critere 
MAP. Cette ambiguite peut etre levee grace a 
1' utilisation des symboles de reference connus du 
recepteur. Cependant, ceci se revele souvent 
insuffisant, car cette densite de probabilite 
conditionnelle possede egalement des maxima locaux 
pouvant etre attaints par I'algorithme SAGE a la place 
de 1' unique maximum global. Pour eviter ce probleme, il 
faut utiliser les symboles de reference pour determiner 

convenablement les conditions initiales {g'^°^},^J . 

L'algorithme SAGE reestime successivement , par 
induction, chacun des vecteurs en gardant inchangees 
les estimations precedentes des autres vecteurs C" , 
avec m^l . Conformement aux proprietes de I'algorithme 
SAGE, ces reestiraations successives garantissent une 
croissance monotone de la densite de probabilite 

conditionnelle a posteriori p({G'},^J!R) - 

Etant donne le vecteur R d' echantillons regus 
durant une tranche temporelle a traiter, le recepteur 
commence par calculer des conditions initiales {g'^^-"},^;, 
des vecteurs ['^'Y/ll ^ partir des echantillons ne 
dependant que des sequences de reference. A cette fin, 
les sequences de reference doivent etre telles qu' une 
separation rapide et simple des contributions des L 
trajets a chacun de ces echantillons soit garantie. 
Pour pouvoir effectuer une estimation initiale, 
permettant le lancement de 1 '/ algorithme d' estimation de 
canal, il faut recourir a une ou plusieurs sequences de 
reference ayant chacune de bonnes proprietes 
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autocorrelation. A I'instar du systems GSM, . ces 
sequences peuvent etre construites a partir de 
sequences CAZAC par extension ou permutation 
circulaire. Ces sequences sont composees de symboles de 
module constant et possedent une fonction 
d' autocorrelation circulaire nulle partout a 
1' exception de I'origine. Dans le cadre de la presente 
invention, il est preferable d' avoir des sequences 
CAZAC avec des symboles appartenant a la constellation 
de la modulation utilisee pour vehiculer les symboles 
de donnees. A titre d'exemple, pour une modulation 
MDP2, la seule sequence CAZAC (a un signe et a une 
permutation circulaire pres) la plus courte, a savoir 
{ -1, +1, +1, +1 } , ne compte que quatre symboles. Pour une 
modulation MDP4, utilisant la constellation , la 

seule sequence CAZAC {toujours a un facteur 
multiplicatif i, -1 ou -i , et a un decalage circulaire 
pres) la plus courte, a savoir , n'est composee que 

de 2 symboles. 

Une sequence de reference basee sur une sequence 
CAZAC de longueur M ne permet d'estimer que M trajets 
du canal discret au plus. Elle est obtenue par une 
extension cyclique de L-1 positions de cette sequence 
CAZAC et comporte au total M+L-1 symboles. La longueur 
M de la sequence CAZAC doit done etre superieure ou 
egale au nombre de trajets L du canal discret a 
estimer. Dans le cadre de la presente invention, il est 
parfois souhaitable d'utiliser des sequences CAZAC 
courtes pour nombre de trajets L a estimer petit 
(inferieur ou egal a 4 pour une modulation MDP2 et a 2 
pour des modulations MDP4, MDP8, ...), car elles 
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permettent de reduire signif icativement le nombre total 
de symboles de reference par tranche temporelle, tout 
en garantissant un echantillonnage plus regulier et 
plus frequent du canal discret, surtout dans les 
5 environneraent a grands etalements Doppler. 

La figure 9 illustre 1' estimation d'un canal 
discret compose de 2 trajets (L=2) dans le cas d' une 
modulation MDP2 . La premiere ligne. A, donne la 
sequence de reference, la seconde, Tl, correspond au 
10 premier trajet, la troisieme, T2, au second trajet et 
la quatrieme, S, donne la composition des echantillons 
regus. La sequence de reference elementaire utilisee 
{-1, +1, +1, +1, -1}, composee de 5 symboles (M+L-l=5) 
n' est autre que 1' extension circulaire d' une seule 
15 position de la sequence CAZAC {-1, +1, +1, +1} composee 
de 4 symboles (M=4). Comme I'indique la figure 5, la 
sequence de reference elementaire occupe les positions 
entre n-4 et n de la tranche temporelle a traiter. Les 
echantillons regus, dependant uniquement de cette 
20 sequence, occupent quant a eux les positions entre n-3 
et n. lis sont donnes explicitement par : 

i?„_3 = Cl,A„,, + Ci,A„_, + W„_, = C°_3 - C^3 + W„_, 
K-i = C2^«-2 + Cl,A„_, + r„_2 = Cl, + C\_, + 
i?„_, = C + C\_^A„_^ + W„_, = Cl_, + C\_, + W„_, 
= C°A„ + C\A„_, +W„= -C° +Cl+W^ 

Pour une estimation peu biaisee des coefficients 
25 C° et Cl du canal discret normalise, pour k pris entre 
n-3 et n, il est necessaire / de faire I'hypothese d'un 
etalement Doppler suffisamment faible pour garantir une 
variation negligeable de canal entre les positions n-3 
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et n . Dans ce cas, une estimation approximative de ces 
coefficients est obtenue en correlant les echantillons 
i?j , k pris entre n-3 et avec la sequence CAZAC 

elle meme et un decalage circulaire d' une position vers 
la gauche de celle-ci. La correlation avec la sequence 
CAZAC {-1, +1, +1, +1} donne : 

C^_3 = Cl_2 = C^_i =Cl= r-i?„_3 + R„_2 + i?„_i + i?J/4 = 

comme premiere estimation des echantillons ii-3,n-2, n-1 
et n du deuxieme trajet^ ou est un terme qui tient 

corapte du bruit et du biais. La correlation avec la 
sequence obtenue apres decalage circulaire, {+1, +1, 
+1, -1}, donne quant a elle 

C°_3 = C„°_2 = Cl, = C° = (R„_^ + R„_2 + K-i - = 
f C^3 + C°_2 + + c° ; / 4 + JV^ 
comme premiere estimation des echantillons aux instants 
n-3,n-2,n-l et n du premier trajet, ou W° est un 
terme qui j oue le meme role que . 

Dans le cas general, 1 ' est imat ion grossiere des L 
coefficients du canal, a un instant donne de la tranche 
temporelle, est obtenue par correlation des 
echantillons dependant uniquement de la sequence de 
reference avec L versions decalees de la sequence CAZAC 
correspondante . 

Chaque sequence de reference construite par 
extension circulaire d' une sequence CAZAC permet 
d'obtenir un echantillon bruite de chacun des trajets 
du canal discret. Le nombre 'd' echantillons par tranche 
temporelle necessaires a la reconstruction du canal 
depend de I'etalement Doppler normalise BdTs et du 
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nombre ISl de symboles composant cette tranche. Plus le 
produit BdTsN est grand, plus le canal varie d'une 
extremite a 1' autre de la tranche temporelle et plus le 
nombre de sequences de references elementaires doit 
5 etre grand. Une sequence de reference elementaire est 
un sous-ensemble de symboles groupes de la sequence de 
reference associee a la tranche temporelle obtenu a 
partir d'une extension circulaire d'une sequence CAZAC. 
Conformement aux proprietes de la decomposition de 

10 KARHUNEN-LOEVE du canal discret, plus une valeur propre 
/I* est grande pour un trajet / donne, plus la 
contribution G^B"' du vecteur propre associe a la 
representation adequate de ce trajet est importante en 
moyenne. La figure 10 montre la repartition des valeurs 

15 propres les plus importantes pour une tranche 
temporelle de N=100 symboles de meme energie normalisee 
a 1 et des etalements Doppler normalises BdTs egaux 
respectivement a 1/50 (courbe 113) , 1/100 (courbe 112) 
et 1/500 (courbe 111) . Cette figure montre clairement 

20 qu'une representation convenable de chacun des trajets 
du canal discret peut s'appuyer, sans perte perceptible 
de performance, sur les quelques vecteurs propres de 
plus grandes valeurs propres. Pour des valeurs faibles 
du produit BdTsN, le vecteur propre de plus grande 

25 valeur propre suffit a lui seul a une representation 
precise de chacun des trajets du canal discret. Ce 
vecteur propre varie alors tres lentement au cours du 
temps et permet done, grace a une sequence de reference 
groupee au milieu de chaqUe tranche temporelle, de 

30 reconstruire un canal discret quasi constant sur toute 
la tranche temporelle. Ceci est le cas du systeme GSM 
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ou les etalements Doppler normalises sont supposes etre 
tres faibles et le canal est suppose etre constant sur 
une tranche temporelle donnee . 

Pour des etalements Doppler normalises un peu plus 
eleves, les deux vecteurs propres associes aux deux 
plus grandes valeurs propres sont necessaires a une 
reconstruction precise du canal discret. Ces deux 
vecteurs varient suffisamment pour qu'une sequence de 
reference composee de deux sequences elementaires , 
placees de part et d' autre du milieu de la tranche 
temporelle, se revele necessaire. 

Dans le cas general, plus grand est le produit 
BdTsN, plus grand est le nombre de vecteurs propres a 
utiliser dans la representation adequate du canal 
discret. De plus, plus grand est ce produit, plus 
importantes sont les variations de ces vecteurs propres 
et done plus important est le nombre de sequences de 
reference elementaires a mettre en oeuvre. De surcroit, 
a longueur de sequence de reference constante, plus le 
nombre de sequences de reference elementaires est 
grand, plus la longueur de chacune doit etre diminuee. 
Pour un nombre de trajets L donne, la longueur d' une 
sequence elementaire, M+L-1, peut au mieux etre reduite 
a 2L-1 en choisissant une sequence CAZAC de longueur 
L ; d'ou I'utilite d' avoir des sequences CAZAC tres 
courtes. Selon 1' invention, la sequence CAZAC {-1, +1, 
+1, +1} de longueur M==4 est recoramandee pour toutes les 
modulations de phase quand le nombre de trajets L est 
inferieur ou egal a 4 et ^ que le produit BdTsN est 
faible. Selon 1' invention egalement, la sequence CAZAC 
{l,i} de longueur M=2 est conseillee pour un nombre de 
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trajets inferieur ou egal a 2 et des produits -BdTsN 
assez importants . 

Soit K le nombre de sequences de reference 
elementaires par tranche temporelle a traiter. Soient 
Cl_^^^,Cl_f^^T...,Cl, / = 0,1,...,I-1 , les echantillons du canal 
discret procures par la A:eme sequence de reference 
elementaire et S,, = {k - M -\-\,k - M -\-2,...,k} I'ensemble des 
indices de ces echantillons. Les conditions initiales, 
{'^"^^^Yi-l ' ^® 1^ representation adequate du canal selon 
1' invention peuvent etre obtenues a partir de ces 
echantillons en utilisant les equations : 



0^=^':^=^ — ou K=im\- et . 

Le facteur de ponderation w[' depend de la norme 
N[ de la restriction des vecteurs propres aux seuls 
echantillons en correspondance directe avec les 
echantillons de canal procures par la sequence de 
reference. II depend egalement de la qualite de la 
contribution du vecteur propre B'* mesuree par le 
rapport A'* / . Plus la variance A* du coefficient G[ 
dans la representation du canal discret est grande par 
rapport a la variance des echantillons de bruit, 

plus la ponderation correspondante est importante et 
s'approche de 1. 

En se basant toujours sur le vecteur R 
d' echantillons regus durant ^une tranche temporelle, le 
recepteur selon 1' invention conduit ensuite une 
estimation iterative du canal discret en partant de ces 
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conditions initiales {g'^°^}^';J . Conf ormement a 
1' extension SAGE de 1 ' algorithme EM, I'estimateur de 
canal optimal salon le critere du MAP calcule 
iterativement, trajet apres trajet, la reestimation 
5 ^ partir de 1' estimation {G'^'^}'^^ . A chaque 

iteration, I'estimateur de canal calcule de proche en 
proche les vecteur G'^''"'^ / = 0,I,...,Z - 1 , a partir des 
vecteurs {g'"'-'^'^}'~1^ et {0"'^"^}'^^^ calcules auparavant. 
Selon 1' invention, la ^eme composante de la (d + l)eme 
10 reestimation g'^^*'^ du vecteur G' associe au /erne trajet 
et obtenue en se servant de 1' expression : 

ou 

et C"*''^ est la « erne composante de la reconstruction, 
C"'^'^^ selon la decomposition de KARHUNEN-LOEVE du meme 
trajet du canal discret a la J eme iteration. Cette 
reconstruction est donnee explicitement par 

2 0 ^^^^mWg™^ _ 

Corame precederament , les ponderations w[ tiennent 
compte de la qualite de la contribution du vecteur B''^ 
a la representation adequate 'du canal discret. Plus la 
contribution d' un vecteur propre est importante par 
25 rapport a la variance du bruit, plus sa ponderation est 
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proche de 1 et plus il joue un role important dan-s la 
reconstruction de 1 ' est imation du canal. 

Afin d' estimer la representation adequate, G'^'^'^'^ 
du /erne trajet en accord avec 1' expression precedente, 
le recepteur commence par calculer une estimation 

C*"''^ A-.- S C*'^ 

de CIA^_, en retranchant les contributions des autres 
trajets m^l . Cette soustraction est faite dans 
I'hypothese ou le vecteur A„ = (^„_(^_,),^„_(^_2),...,^J est 
emis, d'ou I'utilite de ponderer plus tard cette 
expression par la probabilite condit ionnelle a 

posteriori p(p^}^{G"''^'*'^]'1\\G'"^'^Y~'^^ d' avoir reellement 
erais ce vecteur. Une fois cette estimation de C„^„_, 
realisee, le recepteur la multiplie par le conjugue 
de A„_, pour s'affranchir de la modulation. Cette 
operation est effectuee en supposant toujours que le 
vecteur A„ a ete emis, d'ou la necessite de ponderer 
1' estimation resultante : 

de Cl par la probabilite conditionnelle a posteriori 
d' avoir emis ce vecteur. L' estimation finale 

de C'„ est bien sur obtenue en prenant en compte la 
contribution de toutes les valeurs pouvant etre prises 
par A„ . 
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Une fois cette estimation , du « erne echantillon 
du / eme trajet obtenue pour n entre / est N + l-\, il 
ne reste plus au recepteur qu ' a la projeter sur chacun 
des vecteurs B'^ de la base orthonormale associee a la 
representation adequate de ce trajet. La projection sur 
le A: erne vecteur B'^ est ensuite ponderee par w[ afin de 
tenir compte de la qualite de sa contribution a 
1' estimation du /erne trajet du canal a la (<i + l)eme 
iteration. 

En pratique, seuls quelques vecteurs de base 
interviennent de maniere importante dans la 
representation de chacun des trajets. Dans ce cas, le 
nombre de vecteurs utilises dans la transformation de 
KARHUNEN-LOEVE et la transformation inverse peut etre 
reduit sensibleraent sans perte perceptible de 
performances tout en garantissant une reduction 
importante de la complexite. 

L' estimation iterative du canal discret peut etre 
accomplie un nombre fini D de fois de telle maniere 
que 1' estimation globale {c'^°^}^^/J (obtenue a partir de 




par la transformation de KARHUNEN-LOEVE) 



garantisse une degradation imperceptible des 
performances du recepteur par rapport a la solution 
optimale qui serait obtenue apres un nombre infini 
d' iterations . Cette estimation iterative suppose la 
connaissance des probabilites conditionnelles a 

posteriori p(aJr,{g'"^^"'>}^"^'^,{9'"(''^};^;J j pour toutes les 
valeurs possibles des vecteurs A„ = (y^„_(^_,),y4„_^^_2j,...,y4J , 
n = Q,\,...,N + L-2 . Ces probabilites conditionnelles 
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peuvent etre calculees exactement par 1 ' algorithme de 
BAHL, connu aussi sous le nora d' algorithme BCJR [6] . 
Elles peuvent aussi etre calculees de maniere 
approximative et simple grace a des variantes 
linearises de I'algorithme de BAHL [10], 

A titre d' illustration, on presente ici les 
grandes lignes du calcul par I'algorithme de BAHL des 

probabilites conditionnelles p|AjR,{G'"<''*"}|j^g,{G'"^''^}^^J j . A 

cette fin, il faut remarquer, conformement a la figure 
8, que 1 ' echantillon i?„ obtenu a 1' instant n depend a 
la fois du syrabole emis a cet instant et des L-1 

symboles A_p^„-2'---'^«-(l-i) 1® precedant. Ainsi, 
1' ensemble modulateur-canal multi-tra j et s-demodulateur 
peut etre modelise par une chaine de MARKOV a temps 
discret. Cette chaine de MARKOV est decrite par un 
treillis dont les etats representent a chaque instant n 
le contenu du registre a decalage de la figure 8 et 
dont les etiquettes des branches representent la 
version non bruitee des echantillons possibles a la 
sortie du filtre blanchisseur . A titre d'exemple, la 
figure 11 montre une representation en treillis d'un 
canal discret compose de L = 2 trajets pour une 
modulation MDP2 a deux etats {-1, +1} et une longueur 
de tranche temporelle N=4 . 

Un tel treillis represente la sortie attendue du 
canal en fonction du temps et selon les symboles emis. 
Dans le cas illustre ou le canal est suppose presenter 
deux trajets, la sortie du canal est conditionnee par 
le symbole emis a 1 ' instant considere et a 1 ' instant 
anterieur. Par exemple, en partant de 1 ' instant n=l, on 
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peut emettre soit -1, soit +1, sachant qu'a 1 ' instant 
precedent 0, on n ' a rien emis (ncsud initial a 1' extreme 
gauche). Si 1 ' on a emis -1, on s'attend en sortie du 
canal a recevoir (-1) Co + (0)c|) , soit -Cg, quantite qui 
5 figure sur la premiere branche montante vers le nceud 
-1. Si 1 ' on a emis +1, on s'attend a recevoir 
(+1) Cq+0 c'o , soit Cq , quantite qui figure sur la 
premiere branche descendante vers le nceud +1. 

Pour n=l, si 1 ' on emet -1 et si 1 ' on avait emis -1 

10 a 1 ' instant anterieur n=0, on s'attend a recevoir 
(-1) cf+(-l)CJ , soit -C°-c} . Si I'on emet +1, on 
s'attend a recevoir ( + 1) cf-c[=c!'-c[ . De meme, si I'on 
avait emis +1 a 1' instant anterieur, on s'attend a 
recevoir -cf+cl si I'on emet -1, et C°+cJ si I'on emet 

15 +1. Ces quatre valeurs figurent sur les quatre branches 
a n=l. La tranche s ' acheve par une emission de 0, et la 
quantite attendue est -C4 ou C4 selon qu ' on a emis -1 
ou +1 a I'instant anterieur n=3 . Le treillis represente 
done bien, a chaque instant, les quantites attendues en 

20 sortie du canal, selon les symboles envoyes. 

Conformement a une telle representation du canal 
discret, chaque branche du treillis determine une 
valeur unique du vecteur A„ = (^„_(^_i),/i„_(i_2)'---' A) ' car, 
d'une part, les L-1 premieres composantes sont donnees 

25 par I'etat de depart S„ = (^„_(^_,),^^_(i_2),...,-4„_,) de la 
branche et, d' autre part, les L-1 dernieres composantes 
sont donnees par 1/ etat de destination 

S„^j = (4_(i_2)>4-(i-3)'---'^n) • Le probleme de calculer les 
probabilites conditionnelles a posteriori 
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i'|A„|R,{G'"<''^'^}^J^,{G'"^''^}^jj se ramene done au calcul des 
probabilites des branches du treillis 

conditionnellement au vecteur d' echant illons regus R 
et a 1' estimation intermediaire du canal discret 
normalise donnee par les vecteurs c'"'''^" pour les / 
premiers trajets m = 0,l,...J-l et les vecteurs C"^''^ pour 
les L-£ derniers trajets w = /,/ + 1,.. ., Z - 1 . A cette fin, 
les branches du treillis doivent etre etiquetees par 
les coefficients de cette estimation intermediaire et 
non par ceux du canal discret reel. 

L'algorithme de BAHL permet de calculer 
efficacement les probabilites conditionnelles a 
posteriori des branches du treillis en factorisant 
1' expression de la probabilite du vecteur regu R, 
conditionnellement a chacune des sequences A pouvant 
etre emises, afin de ne pas dupliquer toutes les 
operations communes a des sequences emises ayant le 
meme passe ou le meme avenir a un instant donne. A 
cette fin, cet algorithme fait appel a des fonctions 
«„(S) et J3„(S) qui donnent respect ivement un apergu du 
passe et de 1' avenir pour chaque instant n=0, 1, 
N+L-1 et chaque valeur S possible de I'etat S„ . 

La fonction a„(S) est la probabilite conjointe 
d'etre a I'etat S„=S a I'instant n et d'avoir regu les 
echantillons {^,}"J . La factorisation des calculs pour 
des sequences ayant un passe commun permet de calculer 
cette fonction grace a la re9urrence : 

S' 
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ou ;k„(S',S) est la fonction qui donne la probab-ilite , 
conditionnellement a un etat de depart du treillis 
S„.i=S', d'aller a I'etat S„ = S et de recevoir en meme 
temps 1' echantillon R„_^ . Cette expression montre que 
5 1' information sur le passe peut etre acquise a la fois 
a travers les probabilites associees a 1' instant 
precedant immediatement 1' instant actual et a travers 
les probabilites de transitions des branches du 
treillis . 

10 La fonction ^„{S) a une signification equivalente a 

celle de a^{S) . C est la probabilite d' avoir regu les 

echantillons r conditionnellement a I'etat de 

depart S„ = S a 1' instant n. Cette fonction peut etre 
aussi calculee de proche en proche a partir de I'avenir 
15 par la recurrence 

A(S) = ZA.i(S)r„.,(s,s'). 

S' 

La encore, cette recurrence prouve que la 
connaissance de I'avenir pour un etat S„ = S peut etre 
obtenue a partir des caracteristiques a la fois des 

20 etats S^^, =S' suivant immediatement cat etat et des 
branches du treillis les reliant ensemble. Pour une 
branche du treillis caracterisee par I'etat de depart 
S„ = S' = et I'etat d' arrivee 

S„^, = S = (v4„_(.^_,^,...,^„_i, J^) , la fonction ;k^(S',S) est donnee 

25 par : 
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ou -P(^„) est la probabilite a priori d'envoyer le 
symbole A„ . Ces probabilites a priori peuvent 
caracteriser un codage eventual des symboles modules 
reellement emis. lis peuvent egalement caracteriser les 
5 symboles de reference qui sont connus du recepteur et 
qui prennent une seule valeur possible. Pour ces 
symboles de reference, la probabilite a priori -P(^„) 
prend la valeur 1 pour la valeur reellement emise et 0 
pour les autres valeurs possibles de la modulation. 

10 Afin d' ef f ectuer les recurrences sur les fonctions 

a„(S) et J3„(S) f il faut prendre comme valeurs 
initiales : 

ao(So =(0A...,0)) = 1 et ;5^,,_,(S^,,_, =(0A...,0)) = 1 
A partir du calcul par I'algorithme de BAHL des 

15 fonctions a„{S) et y?„(S), il est facile d' obtenir les 
probabilites des branches du treillis, 

conditionnellement aux echantillons regus et a 
1' estimation intermediaire du canal discret, grace a la 
formula : 

oia A = (4_(,_,),4_(,.2),...,^J, S'=(^„_(,_,,,A-(L-2)--.^«-.) et 
S = ),..., ^„_i,^„) • 

II faut remarquer toutefois qua la complaxite da 
I'algorithme de BAHL reste toujours assez elavee at 
25 qu'en pratique il vaut mieux recourir a des algorithmas 
approximatif s qui permettent da diminuar notablemant la 
complaxite des calculs ^ans trop degradar las 
performances. Le recours a ces algorithmes 
approximatif 3 est d' autant plus conseille que le nombre 
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de trajets ou/et d' etats de la modulation . est 
important . 

Apres D iterations, le recepteur dispose d'une 
estimation du canal discret {c'^^^l^^jj proche de la valeur 
optimale selon le critere du maximum a posteriori qui 
est normalement obtenue apres un nombre infini 
d' iterations . A partir de cette estimation du canal et 
des echantillons regus durant la tranche temporelle a 
traiter, le recepteur peut fournir soit des decisions 
fermes grace a 1 ' algorithme de VITERBI, soit des 
decisions ponderees grace a I'algorithme SOVA ou a 
I'algorithme de BAHL. 

L'algorithme de VITERBI cherche la sequence emise la 
plus vraisemblable. II permet done de reduire la 
probabilite d' erreurs de decision sur ces sequences 
emises. L'algorithme SOVA est un algorithme de VITERBI 
avec sorties ponderees. II permet une caracterisation 
plus precise de la fiabilite des decisions effectuees 
par I'algorithme de VITERBI. L'algorithme de BAHL 
permet de determiner la densite de probabilite a 
posteriori de chaque symbole emis, conditionnellement 
aux echantillons regus et a 1' estimation de canal 
obtenue. II permet done de determiner la valeur la plus 
vraisemblable de chaque symbole emis et done de reduire 
le taux d'erreur. Les sorties ponderees sont obtenues a 
travers les probabilites des symboles 

conditionnellement aux echantillons re?us. 

Dans tous les cas, les sorties ponderees ne sont 
utiles que si un decodage ciorrecteur ou un codage de 
source avec entrees ponderees est utilise afin 
d'augmenter les performances du recepteur. 
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A titre d' illustration, dans le cadre. de 
I'algorithme de BAHL, le recepteur utilise 1' estimation 
{C"~°^]^'^ du canal discret pour calculer les fonctions 
a„(S) et P„{S) . Ces fonctions permettent alors de 
5 determiner pour chaque symbole A„ sa distribution de 
probabilite par la formule : 

Toutes les operations precedentes d' estimation 
iterative du canal discret et de detection sont 

10 illustrees sur la figure 12. L' estimateur de canal 120 
utilise des moyens 121 definissant les bases {b'^}^^q , 
I = 0,\,...,L-\ , et des moyens 122 calculant les 
coefficients de ponderation {K'}a=o ®^ WK=o " 
L' estimateur 120 met en oeuvre une premiere estimation 

15 grossiere du canal a partir des seules sequences de 
reference elementaires par correlation circulaire avec 
les sequences CAZAC correspondantes . Cette premiere 
estimation est symbolisee par le bloc reference 
1201, qui delivre la representation adequate de 

20 1' estimation initiale • Ensuite, 1' estimateur 

120 met en oeuvre une reestimation d'ordre <i , 
representee par le bloc 3^ reference 1202, qui delivre 
la representation adequate de 1' estimation de canal 
{g'*'^^}^^_"J , suivie d' une reestimation d' ordre d + lr 

25 representee par le bloc 5^^, reference 1203, qui delivre 
la representation de 1' estimation de canal {g^''^"'"}^^^ . 
Enfin, 1' estimateur 120 met en oeuvre une derniere 
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reestimation d'ordre D, representee par le bloc 3^ 
reference 1204, qui delivre la representation adequate 
de 1' estimation de canal {g"'°^}'^~^ . A partir de cette 
representation adequate du canal estime, I'estimateur 
5 fournit par les moyens 1205 1' estimation definitive 
I ^^j-J canal discret en raettant en oeuvre la 

transformation de KARHUNEN-LOEVE . 

La fourniture a un decodeur correcteur d' erreur ou 
a un codeur de source de decisions fermes ou souples 
10 sur les symboles emis, s'effectue par des moyens 123 
qui mettent eventuellement en oeuvre 1' algorithme de 
VITERBI [3], 1' algorithme SOVA [2], 1' algorithme de 
BAHL [6] ou I'algorithme de BAHL linearise [10]. 

Selon 1' invention, la reestimation iterative du 
15 canal est effectuee D + 1 fois et permet de calculer 

{g'^^^ITo'' {G'^'^I^TJ. {G'^''^}/to' successivement. La 

figure 13 illustre les L etapes composant chaque 
iteration et qui sont necessaires a 1' estimation 
successive de chacun des L trajets. Grace aux moyens 

20 121 et 122 de la figure 12, (fournissant les bases 
{B^^'I^Jq et les ponderations {^"liJ ) ' (<i + l)eme 

iteration 3d+i de I'estimateur de canal met en oeuvre 
une premiere etape, basee sur 1' algorithme de BAHL et 
representee par le bloc BAHLq reference 130, pour 

25 delivrer les probabilites condit ionnelles a posteriori 
p|A„|R,{G"'^''^}^^'|jj a partir de la d erne reestimation du 
canal. A partir de ces probabilites conditionnelles , la 
(cf + l)eme iteration met ensuite en oeuvre une deuxieme 
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le bloc SAGEq reference 131, pour permettre la ( J + 1) erne 

reestimation G°'''^'^ du premier trajet. 

A partir de la (<i + 1) erne reestimation des I -I 
premiers trajets, la (J + 1) erne iteration met en oeuvre 
une troisieme etape, basee de nouveau sur I'algorithme 
de BAHL et representee par le bloc BAHL, reference 132, 
pour delivrer les probabilites conditionnelles a 

posteriori p[AjRlG'"^'''^}'2A^'"^'^}l'',) ^ Partir de la d erne 
reestimation des L-l derniers trajets du canal et de 
la (cf + i)eme reestimation des / premiers trajets de ce 
canal. Par la suite, la (J + l)eme iteration met en 
oeuvre une quatrieme etape, basee sur I'algorithme SAGE 
et representee par le bloc SAGE/ reference 133, pour 
permettre la (cf + l)eme reestimation g'^''^" du /erne 
tra j et . 

Enfin, a partir de la (d + l)eme reestimation des 
L-l premiers trajets, la {d + \)eme iteration met en 
oeuvre une derniere etape, basee sur I'algorithme de 
BAHL et representee par le bloc BAHL^_, reference 134, 
pour delivrer les probabilites conditionnelles a 
posteriori p(a„|r,{G'"^''"'^}^;|,,G^'-'^'<''' j a partir de la £^ erne 
reestimation du dernier trajet du canal et de la 
(cf + l)eme reestimation des L-\ premiers trajets de ce 
canal. Par la suite, la (d + \)eme iteration met en 
oeuvre une derniere etape, basee sur I'algorithme SAGE 
et representee par le bloc SAGE^., reference 135, pour 
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permettre la (cf + 1) erne reestimat ion QU-m^^^) -leme 
tra j et . 

La premiere iteration, symbolisee par le bloc 3o 
reference 1201 de la figure 12, ne peut pas mettre en 
oeuvre les algorithmes de BAHL et SAGE pour calculer la 
premiere estimation du canal discret. Selon 
1' invention, on utilise alors les symboles de reference 
dont les probabilites conditionnelles a priori et a 
posteriori sont bien connues du recepteur pour 
effectuer cette premiere estimation. Cette premiere 
estimation s'effectue, en correlant les echantillons 
regus, dependant uniquement de ces symboles de 
reference, avec les sequences CAZAC correspondantes . 
Elle met en oeuvre les moyens 122 de la figure 12 pour 

calculer les ponderations • 
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REVENDICATIONS 

1. Procede de communications radiomobiles du type 
a acces multiple a repartition en temps, dans lequel on 
alloue a plusieurs utilisateurs des tranches 
temporelles pour 1' emission et la reception de signaux 
radioelectriques, et dans lequel : 

- a 1' emission, on constitue, dans chaque tranche 
temporelle allouee, des symboles de donnees et 
des symboles de reference, on emet un signal 
radioelectrique contenant ces symboles de 
donnees et ces symboles de reference, ce signal 
empruntant un canal radiomobile mult i-tra j ets , 

- a la reception, on effectue un filtrage adapte 
du signal regu, on echantillonne le signal 
filtre au rythme des symboles, on effectue une 
estimation du canal radiomobile, on traite les 
echantillons en tenant compte de cette 
estimation du canal et 1 ' on restitue les 
symboles de donnees propres a la tranche 
allouee, 

ce procede etant caracterise en ce que : 

- a 1' emission, pour chaque tranche temporelle 
allouee, on constitue plusieurs sequences de 
reference possedant des proprietes de correlation 
temporelle appropriees et 1 ' on repartit ces 
sequences dans chaque tranche temporelle, 

— a la reception : 

i) pour estimer le canal radiomobile, on 
represente ce canal par un treillis 
comprenant des branches symbolisant les 
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transitions entre deux etats successifs du 
canal et caracterisant les sequences 
possibles i la sortie du canal, 
ii) on calcule les probabilites de ces branches 
5 en mettant en ceuvre un algorithms -iteratif 

dans lequel on calcule, par un nombre fini 
(D+1) d' iterations et selon le critere de 
probability maximum a posteriori, une 

estimation d-' une representation des 
10 L trajets de ce canal, 1' iteration de rang 

(d+1) (3d+l) permettant d'obtenir une 

reestimation {^'^''^'^K ^ partir de 
1 ' estimation (g^^'^^I^Jo obtenue par 1 ' iteration 
precedents {3d), et, a chaque iteration 

15 (5d+l) , on calcule en L etapes les vecteurs 

Qe(d+i)^ ^=0, 1, L-1 a partir des 

vecteurs ^"'^'^^^}l=o et ^"^^'^^1, calcules a 
I'etape precedente, 1 ' estimation |c'^'^^}fj des 
L trajets du canal etant final ement obtenue 

20 a partir de 1' estimation de sa 

representation ^^^^^ft^l . 

2. Precede selon la revendication 1, dans lequel, 
a chaque etape de chaque iteration, on met en oBuvre 
25 successivement : 

i)un algorithme de BAHL ou une de ses variantes 
pour deiivrer des probabilites 

conditionnelles a partir de la d^^"*® 
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reestimation des L-£ derniers trajets et de 
la (d+l)^^™* reestimation des i premiers 
trajets, 

ii) un algorithms SAGE pour delivrer une (d+l)"-^^ 
5 reestimation du i^^'^^ trajet. 

3. ProG^de selon la revendication 1, dans lequel, 
ayant obtenu une estimation du canal discret |c'^^*};^o 
et, compte tenu des echantillons regus durant la 

10 tranche temporelle traitee, on prend une decision yur 
la sequence des symboles emis en utilisant un 
algorithme optimal de BAHL, ou une de ses variantes, 
pour obtenir la probabilite a posteriori de chaque 
syinbole emis. 

15 

4. ProcedS selon la revendication 1, dans leqiael 
on r^partit K sequences de reference dans chaque 
tranche allouee avec K compris entre 2 (inclus) et 12 
(inclus) . 

20 

5. Procede selon la revendication 1, dans lequel 
chaque sequence comprend un nombre de symboles de 
reference inferieur k 8. 



S. Procede selon I'une quelconque des 
revendications 1 a 5, dans lequel les sequences de 
reference possedent une fonction d ' autocorrelation 
circulaire nulle partout sauf a I'origine. 

7. Procede selon la revendication 6, dans lequel 
on utilise une sequence de reference {-1, +1, +lr +1} 
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ou une sequence derivee de celle-ci par pe2miut;ation 
circulaire, ou changement de signe pour un nombre de 
trajets (L) inferieur ou egal a 4. 

8, proced^ selon la revendication 6, dans lequel 
on utilise une sequence de reference {1/i} ou une 
sequence derivee de celle-ci par permutation 
circulaire, conjuguaison ou changement de signe, pour 
un nombrs de trajets (L) inferieur ou ^gal a 2. 
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